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Еволющйна модель задач! булева програмування 


В работе представлены результаты исследования фрагментарной и эволюционной модели задачи булева 
программирования. Показано, что при определенных условиях задача булева программирования может 
рассматриваться как задача на фрагментарной структуре. Предложена эволюционно-фрагментарная модель 
задачи на множестве перестановок с геометрическим оператором кроссовера. Метод протестирован на 
наборе индивидуальных задач различных размерностей. 

Ключевые слова: задача о рюкзаке, фрагментарная структура, эволюционная модель. 
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У робот! представлен! результати досллдження фрагментарно! 1 еволющ но! модел! задач! булева 
програмування. Показано, що при певних умовах задача булева програмування може розглядатися як задача 
на фрагментарнй структур!. Запропоновано еволющйно-фрагментарна модель задач! на множин! 
перестановок з геометричним оператором кросоверу. Метод протестований на набор! 1ндив!дуальних задач 
рвних розмрностей. 

Ключов! слова: задача про рюкзак, фрагментарна структура, еволющйна модель. 


Введение 


Во многих прикладных оптимизационных задачах отыскание точного решения 
является достаточно трудной проблемой. Иногда единственным методом, который 
гарантирует желаемый результат, является лишь полный перебор всех допустимых ре- 
шений. Аналогичная проблема возникает при отыскании точного решения многих за- 
дач распознавания, в частности при отыскании решения задач, для которых доказано 
свойство №Р-полноты. 

Для поиска решения подобных задач, как правило, используются приближенные 
и эвристические алгоритмы. Одним из эвристических подходов является подход, ос- 
нованный на применении фрагментарных алгоритмов [1], [2]. Однако использовать 
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фрагментарный алгоритм для поиска оптимального решения задач с четко определенным 
формальным критерием вряд ли разумно. Фрагментарный алгоритм останавливается 
при отыскании первого максимального по включению фрагмента и, следовательно, не 
гарантирует качества решения с точки зрения оптимума критерия. Для решения диск- 
ретных оптимизационных задач оправданным является метод, основанный на комби- 
нации эволюционного и фрагментарного алгоритмов. В настоящей работе подобная 
технология будет применена к известной задаче булева программирования — задаче о 
рюкзаке [3]. 

Постановка задачи. Задачу о рюкзаке будем рассматривать в следующей поста- 
новке: на множестве и-мерных двоичных векторов 


25 =4{х|х=(х,х,,...х,), х, Е {0,1 =1,2,...,пт} 
найти вектор х =(х,,х,,...,х,) , для которого значение линейной целевой функции 
Е(х) = ах +сх, +...+с,х, (1) 
максимально при условии выполнения т линейных неравенств 


ах ал фортах еь 


ета ото та хе в 6) 


2х +а,2х. +...+атх, ЗВ, 


т 


Коэффициенты с,а — неотрицательные вещественные числа, 1=1,2,...,И, 


7 Л 
2 

В [4] доказано, что задача о ранце является МР-полной. Одним из методов 
отыскания точного оптимального решения этой задачи является метод ветвей и 
границ, описание которого можно найти во многих источниках [3], [4]. 

Простой, широко известный жадный алгоритм отыскания приближенного реше- 
ния задачи о рюкзаке заключается в следующем: для каждой переменной х», {=1,2,...,И 
вычисляем ее удельную ценность 


Б, 
у= шш 4{6—} 
аз>0,./=1,2,..т а 
7 
и упорядочиваем переменные х; в порядке убывания удельной ценности ХХХ: 


п 


На начальном шаге полагаем х° = (0,0.....,0). На каждом последующем шаге оп- 
ределяется значение координаты й#, в решении. Изменим решение х', положив коор- 


динату х, равной 1. Если измененное решение удовлетворяет всем ограничениям (2), 


то полагаем у. равным этому измененному решению. В противном случае оставляем 


К К-1 
Хх =х . 


Конечно, приведенный алгоритм не гарантирует точности полученного с его по- 
мощью решения, а является лишь простой эвристикой, которая позволяет найти неко- 
торое допустимое решение задачи. Попытаемся обобщить эту эвристику с помощью 
понятия фрагментарной структуры. 
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Фрагментарная структура 


Определение 1. Фрагментарной структурой (Х,Е) на конечном множестве Х на- 
зывается семейство его подмножеств Ё={ЁЕ,Е,,....ЁЕ,}, такое, что УЕ ЕЁ, Е =03е 


ЕЁ, Е, \{е} ЕЕ. 

Элементы из множества Е будем называть допустимыми фрагментами. Таким 
образом, для любого допустимого фрагмента Ё; существует нумерация его элементов 
Е, ={ел»ер,...е„ }, такая, что У =1,2....,5, {е1,ер,...е,}еЁЕ. 

Определение 2. Одноэлементные множества, которые являются допустимыми 
фрагментами, будем называть элементарными фрагментами. 

Определение 3. Фрагмент называется максимальным, если он не является под- 
множеством никакого другого фрагмента. 

Очевидны следующие свойства фрагментов: 

Свойство 1. Пустое множество является фрагментом, О ЕЁ. 


Свойство 2. Пусть тах |, | = М . Тогда в Ё найдутся элементы, мощности которых 
ТЕГИ 


будут М, М -ЬМ -2....,0. 
Теорема 1. Если (Х, Е) — фрагментарная структура на множестве Х, то для любого 


непустого множества А е Е существует нумерация его элементов А = О ы х, } та- 
кая, что УК, К =1и множество {х,х,,...х, |} Е Е. 

Доказательство. Пусть 4 = А, АЕЁи || =и> 0. Выберем элемент х, е 4, таким 
образом, чтобы 4 = 4 \{х,} ЕР. Повторив эту процедуру и — 1 раз по отношению ко 


множествам 4, 4,,..., А, , получим требуемую нумерацию элементов множества А. 

Из теоремы вытекает, что всякий допустимый фрагмент можно построить из 
пустого множества, последовательно добавляя к нему элементы так, чтобы на каждом 
шаге такой процедуры полученное подмножество было допустимым фрагментом. 

Максимальный фрагмент может быть построен с помощью следующего «жадного» 
алгоритма: 

а) элементы множества Х линейно упорядочиваются; 

6) на начальном шаге выбирается пустое множество Х’, = ©; 


в) на шаге с номером А+1 выбирается первый по порядку элемент хе Х \ Х,,, такой, 
что ХМ еЕ; 

г) алгоритм заканчивает работу, если на очередном шаге не удалось найти элемент 
хе Х \ Х, стребуемым свойством. 


Приведенный выше алгоритм построения максимального фрагмента во фрагмен- 
тарной структуре будем называть фрагментарным алгоритмом. Результат применения 
фрагментарного алгоритма определяется заданным линейным порядком на множест- 
ве Х. Таким образом, любой максимальный фрагмент может быть описан некоторой 
перестановкой элементов множества Х. Пусть АеЁ. Условие для элемента хеХ, 


при котором А ®) {х} ЕЁ ь будем называть условием присоединения элемента х. 


Теорема 2. Если АЕЁ и Ухе Х существует алгоритм полиномиальной трудо- 
емкости по числу элементов множества Х для проверки условия присоединения, то 
задача построения максимального фрагмента является полиномиально разрешимой. 


«Штучний 1нтелект» 2013 №1 125 


4к Козин И.В. 


Доказательство. Пусть трудоемкость проверки условия присоединения состав- 
ляет 0("), где к натуральное число. Трудоемкость упорядочения элементов может 


быть оценена величиной О(иши). Тогда для фрагментарного алгоритма построения макси- 
мального фрагмента трудоемкость оценивается величиной О(пти иен ) ‚ то есть О("" ). 


Покажем теперь, что задача булева программирования в постановке (1) - (2) 
является задачей на фрагментарной структуре. Допустимым фрагментом в этой струк- 
туре является такой набор индексов [< Х = {1,2...., и}, для которого 


о 


1=/ 


Очевидно, что набору индексов / соответствует допустимое решение задачи (1) -— (2) 


1) 1Е!Г | 
х=(х,х.,...Х,) где х, = : 1=1,2,...И. 
0, 11 


Значение целевой функции на этом решении Е(х) = уе . Это же значение будем 
1Е1 


рассматривать как вес соответствующего набора индексов. 
Таким образом, задача булева программирования может быть сформулирована, как 
задача на фрагментарной структуре (Х,{1}) : найти допустимый фрагмент максимального 


веса. Очевидно, оптимальное решение этой задачи можно искать среди максимальных фраг- 
ментов. Если У1=1,2,....п с, >0, то любое оптимальное решение задачи может быть 


получено фрагментарным алгоритмом при соответствующей перестановке элементарных 
фрагментов. Однако число перестановок фрагментов в модели есть и! и потому задача 
перебора всех фрагментов является трудной. Для быстрого поиска решений, которые 
близки к оптимальному, предлагается следующая эволюционная модель на фрагментар- 


ной структуре. 


Эволюционная модель 


Эволюционные (генетические) алгоритмы подробно рассматривались в много- 
численных публикациях [5], [6]. Для ряда оптимизационных задач удалось предложить 
достаточно эффективные процедуры поиска оптимальных решений, основанные на при- 
менении эволюционных алгоритмов. Для реализации эволюционного алгоритма необхо- 
димо выделить ряд объектов и процедур, совокупность которых будем называть эволю- 
ционной моделью [7]. Основные составляющие эволюционной модели следующие: 

— базовое множество решений — множество допустимых решений Х, на котором 
ищется оптимальное решение задачи; 

— оператор построения начальной популяции: процедура, которая позволяет выде- 
лить на множестве всех допустимых решений его подмножество У < Х для последующей 
ЭВОЛЮЦИИ; 

— критерий селекции — алгоритм, который позволяет сравнивать по качеству реше- 
ния в рамках заданной популяции; 

— оператор кроссовера К: ХхХ -> Х ‚, позволяющий по двум допустимым реше- 
ниям-родителям построить новое решение-потомок из множества допустимых решений; 

— оператор мутации М: Х > Х; 

— оператор отбора, который выделяет множество пар в У для выполнения 
операции кроссовера; 
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— оператор эволюции, позволяющий строить новые популяции из множества ро- 
дителей и потомков; 

— правило остановки, которое определяет условие остановки эволюционного алго- 
ритма. 


Отобранные 
пары - 
родители 


Текущая Потомки 


популяция 


Базовое 
множество 


Промежуточная 
популяция 


1. Оператор построения 
начальной популяции 
Оператор селекции 

Оператор кроссовера 
Оператор мутации Выбор Условие 
Оператор эволюции решения остановки 


> 


Рисунок 1 — Эволюционная модель 


Опишем кратко принцип работы эволюционного алгоритма. На начальном шаге с 
помощью оператора начальной популяции строится множество решений У%. На каждом 
очередном шаге предполагается заданным некоторое множество допустимых решений — 
текущая популяция. На первом шаге это множество У =. Для каждого из элементов 


множества У вычисляется значение критерия селекции. Далее с помощью оператора 
отбора в текущей популяции У выбирается множество пар для кроссовера. К каждой 
паре из выбранного множества пар применяется оператор кроссовера, а затем к ре- 
зультату кроссовера применяется оператор мутации. Таким путем находится множество 
элементов — потомков У. 


К промежуточной популяции У у. которая является объединением текущей по- 
пуляции и множества потомков, применяется оператор эволюции, который выделяет на 
этом множестве новую текущую популяцию. Процесс эволюции повторяется до тех 
пор, пока не будет выполнено условие остановки эволюционного алгоритма. Блок- 
схема эволюционного алгоритма приведена на рис. 1. 

Свойства фрагментарных структур позволяют построить особый класс эволюцион- 
ных алгоритмов на фрагментарных структурах — ЭВФ-алгоритмы. 

ЭВФ-алгоритм является комбинацией эволюционного и фрагментарного алгоритма. 
Опишем эволюционную модель и принцип действия такого алгоритма. 

В качестве множества допустимых решений рассматривается подмножество макси- 
мальных фрагментов на заданной фрагментарной структуре. Каждый фрагмент из этого 
множества определяется как результат работы фрагментарного алгоритма при некоторой 
заданной перестановке элементарных фрагментов. Таким образом, любому допустимому 
решению соответствует определенная перестановка чисел 1,2...., М, где № — количество 
элементарных фрагментов. Для каждого допустимого решения определено значение це- 
левой функции. 
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Базовое множество Х эволюционной модели — это множество 5, = {й,5,,..., 1, } всех 


перестановок чисел 1,2,...„№. Оператор построения начальной популяции выделяет под- 
множество заданной мощности О) из множества Х. 

Правило вычисления критерия селекции устроено следующим образом: по за- 
данной перестановке фрагментов с помощью фрагментарного алгоритма строится 
максимальный допустимый фрагмент. Вычисляется значение целевой функции зада- 
чи для этого фрагмента. 

Опишем теперь оператор кроссовера. Пусть И = (и, и,,....иу)и Г =(%,,,....Уу) — 


две произвольные перестановки. Перестановка-потомок строится следующим образом: 
последовательности И и Г просматриваются с начала. На А-м шаге выбирается наи- 
меньший из первых элементов последовательностей и добавляется в новую переста- 
новку-потомок. Затем этот элемент удаляется из двух последовательностей-родителей. 
Например, 


К((2,3,4,7,8,1,6,5),(3,4,6,2,1,5,8,7)) = (2,3,4,6,1,5,7,8). 


Оператор мутации М выполняет случайную транспозицию (замену местами двух 
элементов) в перестановке. 

Оператор селекции выбирает случайным образом набор пар из заданного числа пар 
во множестве перестановок текущей популяции. 

Оператор эволюции элементы промежуточной популяции упорядочивает в после- 
довательность по убыванию значения критерия селекции. В качестве новой текущей 
популяции выбираются первые О элементов последовательности. 

Обычное правило остановки — количество поколений достигло предельной границы /. 
Лучшая по значению критерия селекции перестановка из последней построенной по- 
пуляции определяет приближенное решение задачи. 

Предложенный подход является универсальным и позволяет применять один и тот 
же эволюционный алгоритм к любым оптимизационным задачам на конечных фрагмен- 
тарных структурах. 


Результаты тестирования 


Разумным представляется сравнение работы ЭВФ-алгоритма с другими известными 
алгоритмами. Конечно, хотелось бы получить численные оценки сходимости алгоритма, 
сравнение по скорости и точности получаемого решения. Однако, как правило, для \№Р- 
трудных задач, в которых применяется ЭВФ-алгоритм, это не представляется возможным. 

Описание эксперимента. 

Входными параметрами при описании серии случайных задач являются: число пе- 
ременных и; число ограничений т; диапазон изменения коэффициентов [41,4] в целевой 
функции и ограничениях задачи; количество задач в серии 5. 

С помощью генератора случайных чисел генерируется набор коэффициентов с» @,6; 
в условии задачи о рюкзаке. 

Рассматривались 4 серии задач. Серия А — задачи малой размерности с числом пе- 
ременных и=10, и с числом ограничений 7и=6. Серия Б — с числом переменных 20 и ко- 
личеством ограничений 12. Серия В — задачи с числом переменных 50 и количеством 
ограничений 30. Серия Г — задачи с числом переменных 50 и количеством ограничений 30. 
В каждой серии генерировалось 100 задач. 

Задачи решались с помощью известного приближенного алгоритма — упорядочения 
по максимальной удельной ценности, методом случайного поиска и ЭВФ-алгоритмом. 
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Сравнение алгоритмов осуществлялось по следующим направлениям: 

Рекорд — количество задач в серии, где алгоритм оказывался лучшим среди 
тестируемых. Рейтинг по Борда — сумма числа баллов, набранных на каждой задаче 
серии. За первое место в сравнении назначалось 5 баллов, за второе 4, за третье 3. 

Результаты сравнения алгоритмов на сериях задач приведены в нижеследующих 
таблицах. 

Таблица 1 — Результаты применения различных подходов 


Прибл.алгоритм Случ.поиск ЭВФ-алгоритм 
Серия Число задач 

Рекорды | Рейтинг | Рекорды | Рейтинг | Рекорды | Рейтинг 
А 10х6 100 44 277 100 400 100 400 
Б 20х12 100 10 221 54 348 100 400 
В 50х30 100 И 280 0 204 91 391 
Г 100х60 100 48 347 0 201 54 354 

Выводы 


Теоретические результаты и результаты численных экспериментов показывают, что 
ЭВФ-алгоритм может достаточно эффективно использоваться как эвристический алгоритм 
в при решении оптимизационных задач линейного булева программирования. ЭВФ-алго- 
ритм является управляемым и качество его работы может возрастать при изменении ряда 
параметров алгоритма, таких как, величина популяции, число пар для селекции, количество 
поколений, количество эволюций и т.д. 
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ВЕЗОМЕ 


1.7. Кот 
Еуоинпопагу Моае[5 Воо[еап Ргоэтаттте РгоМет 

\е за4у меП-Кпо\т ргоет оЁ Водеап ргоэгапиише - Кпарзаск рго ет. ТЬ1$ 
ргоет Ваз патегои$ аррИсаНоп$ ш епошеение ап4 есопопсз. 

То Нпа 1е арргохипае зо[айопз$ оЁ Фе Кпарзаск ргоет 15 ргорозе4 ю зе Ше 
еуо[аНопагу-Наотетагу то]. 

Те рарег 4езстоез ш деа! Фе сопсер{ оРа Назтещеа згис@хге. Зоте ргорегйез 
ОГ Настешагу згис@аге$ аге ргоуе4. Егастещагу а1еотИбт, {ай аПо\уз Юг ап ог4ег оЁ 
е]ететагу Настеп 10 сопзбасё а тахипа| аПоууаМе Наотепё оГ а Наотещагу 
угистоге, 15 дезсге4. 

ЗВо\т ай е Кпарзаск рго ет сап Бе зееп аз а ргоМет Ул Фе Настещагу 
угисфиге. То зеагсН г йе арргохипае орНта| зоайоп о 11$ ргоШет 15 ргорозеа 
еуо[аНопагу а!еотИбт, г УТиеЬ Фе Базе ей 1$ Фе зе{ оР ретпщаНопз оЁ @“етещагу 
Каотеп. Тре питенса| ехрегитепе оп а 1агое питЬег оЁ {е5ё ргоМепз зВо\е4 Фе 
еНесНуепе$$ оЁ 1е еуо[аНопагу тоде! Гог а зеагсН оЁ Ве орйта| зо[аНоп оЁ Ве рго ет 


ог Водеап ргоэтатит!е. 


Статья поступила в редакцию 02.11.2012. 
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